QCD EVOLUTION OF HADRON AND JET 
MULTIPLICITY MOMENTS * 



WOLFGANG OCHS 
Max-Planck-Institut fur Physik, Fohringer Ring 6 
D-80805 Miinchen, Germany 
E-mail: wwo@mppmu.mpg.de 

We describe recent applications of the MLLA evolution equations to the calcu- 
lation of mean multiplicities in quark and gluon jets and the higher moments of 
hadron and sub-jet multiplicity distributions in e + e~-annihilation as function of 
c.m.s. energy Q and resolution parameter y C ut- The transition from jets to hadrons 
with increasing jet resolution is considered. 



1 Introduction 

Multiparticle production in hard collision processes is described within QCD by combin- 
ing the perturbative approach to the parton cascade evolution with a non-perturbative 
treatment of the transition towards the hadronic final state. The perturbative phase is es- 
sentially determined by the QCD scale parameter A and, possibly, the quark masses. The 
gluon bremsstrahlung which dominates the partonic cascade process with its collinear and 
soft singularities leads to the characteristic jet structure of the multiparton final state. 
The formation of partonic jets can be quantitatively studied by constructing jets explic- 
itly using an algorithm which combines partons into jets at an externally given resolution 
scale (parameter Q c ). 

In the calculation of jet production phenomena one relates hadronic jets with partons 
at the same resolution scale neglecting in general the effects of hadronization. ft is an 
interesting question down to which scale one can follow this scheme of identifying a 
parton and a hadron jet which we will adress below. 

In a particularly simple ansatz for hadronization one assumes that observables for the 
multiparticle final state are proportional to the corresponding quantities for partons at a 
characteristic resolution scale Qq 1 an idea which has been proposed originally for single 
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inclusive energy spectra and has been called "Local Parton Hadron Duality" (LPHD 
Subsequently, this idea has been applied to a wider range of phenomena including inclu- 
sive multiparticlc correlations and even quasi-exclusive processes (for reviews, see [21 E])- 
The agreement with data generally increases with the accuracy of the calculation. The 
question arises whether this kind of agreement can be derived from a general principle 
or should the agreement be considered as largely accidental. Whereas there is no gen- 
erally accepted answer it is worth while to check this correspondence for more complex 
observables and at the same time to increase the accuracy of the predictions aiming at a 
reliable phenomenology. 

Here we will consider in particular the moments of multiplicity distributions of hadrons 
and jets. We investigate to what extent the jet observables can be connected with the 
corresponding hadron observables in the limit 

Qc^Qa- (1) 

Such a connection would suggest a one-to-one correspondence between hadrons and par- 
tons. Clearly, such a correspondence cannot exist in general for any exclusive limit but 
only after a certain averaging in the sense of a dual description. 

The mean multiplicities and multiplicity moments in quark and gluon jets as function 
of primary energy and (sub-) jet resolution has been calculated using the evolution equa- 
tion in the Modified Leading Logarithmic Approximation (MLLA) This equation 
has been presented originally at the Leningrad Winterschool 1984 [Jj], is explained in 0, 
and can be considered as an extension of the well known DGLAP evolution equation 
towards small particle energies taking into account soft gluon coherence as realized in 
a probabilistic way by angular ordering [7j. Whereas there have been various levels of 
analytical approximations to the solution of the MLLA evolution equation with increas- 
ing number of subleading logarithmic terms we find that the numeric solution of this 
equation yields quantitative agreement with the variety of global observables discussed 
here. 

2 Moments of multiplicity distribution 

Be P n the distribution of particle (parton) multiplicity in a jet. Then we consider the 
unnormalized and normalized factorial moments f q and F q 

oo 

/ 9 = ^n(n-l)...(n-g + l)P„, F q = f q /W, N = fx (2) 

n=0 
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with mean multiplicity N. Furthermore, one introduces the cumulant moments k q and 
K q which are used to measure the genuine correlations without uncorrelated background 
in a multiparticle sample 



9-1 



= 1 



k 1 = h ~ 2^ ■ K-ifh K q = k q /N\ (3) 



in particular K2 = F2 — 1, K3 = F3 — 3i<2 + 2; for a Poisson distribution K\ = 1, K q = 
for q > 1 . The moments can be conveniently computed with the help of the generating 
function 

00 

Z(u) = Y. P ^ ul% W 

n=l 



_ d n Z{u) { _ 9"lnZ( M ) , 



Of special interest are the ratios 

JJ, = /if 9 /F, (6) 

which have been predicted to show an oscillatory behaviour at high energies [S] with 
the first minimum near q ss 5 at the mass of the Z boson. Such a minimum has been 
observed indeed in e + e~ annihilations at SLC (HJ and at LEP ^OIE] but the magnitudes 
of moments are found much smaller than originally expected in |S] . 

3 Perturbative QCD predictions 

Predictions on the global quantities as above can be obtained from the MLLA evolution 
equation for the generation function Z(Y C , u) and the initial condition at threshold which 
read in the simplified world of gluodynamics without quarks 



— z(y c ,u) = r za dz ^L 



{Z{Y C + In z, u)Z(Y c + ln(l -z),u)- Z{Y Cl u)} (7) 
Z(0, u) = u. (8) 

In the general case there are two coupled equations for Z g and Z q . In Q the evolution 
variable is Y c = In M- in the jet energy E and resolution parameter Q c . The parameter Q c 
limits the transverse momentum from below kx = min(z, 1 — z)E > Q c . This restriction 
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yields a parton cascade with a minimum kx separation and can be compared with the 
jet ensemble constructed from hadrons using the so-called kx- or "Durham" -algorithm. 
The initial condition JSJ sets the multiplicity to N = 1 at threshold E — Q c and F q = 
for q > 1. 

Asymptotic solutions can be obtained in the Double Logarithmic Approximation 
(DLA) which includes only the dominant contributions from the collinear and soft sin- 
gularities, i.e. the splitting function P gg {z) ~ 1/z in Q; the next to leading single log 
terms are included in the MLLA. Up to this order the results are complete; further log- 
arithmic contributions beyond NLLO can be calculated, but they are not complete and 
neglect in particular process dependent large angle emissions. Nevertheless they improve 
the results considerably as they take into account energy conservation with increasing 
accuracy. The full solution of Eq. J7J) , corresponding to the summation of all logarithmic 
orders, can be obtained numerically. Alternatively, one may calculate results of the QCD 
cascade from a Monte Carlo generator, we compare here especially with ARIADNE |12j . 
which is based on similar construction principles as Eq. J7J). 

4 Mean particle multiplicity in quark and gluon jets 

The multiplicities N g , N q in gluon and quark jets can be obtained from the MLLA evo- 
lution equations. At high energies one can write 



iv 9 (y)~exp^y y 7 (y)^j 



(9) 



where the anomalous dimension 7 has the expansion in 70 = y/2C \a. s ji: 

7 = 70 (1 - 0170 - a 2 7o - a37o • ■ •): ( 10 ) 
likewise the ratio of gluon and quark jet multiplicity 

JV Ca 

rEE ~N~ = cv( 1 ~ mo ~ r:n o _r 37o---)- (11) 

with the colour factors Ca = 3 and Cf = |. The coefficients Oj and Ti can be obtained 
from the evolution equations. 

The rise of parton multiplicity in a quark jet is given in MLLA by 
N oc exp [ciy/ln(E/A) + C2 lnln(i?/A)] and this formula describes well the data in e + e~ 
annihilation at LEP1 and LEP2 (for reviews, see p3E] and [T5]L 
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Figure 1: The ratio of the mean multiplicities in gluon jets and quark jets N g and N q . 
Results from evolution equations of different order of approximation in comparison with 
experimental data obtained in e + e _ -annihilation. 



The role of higher logarithmic orders can be studied in the behaviour of the mul- 
tiplicity ratio r in The asymptotic limit r = Ca/Cf acquires large finite energy 
corrections in NLLO QUE! and 2NLLO order [El El 

- - 2 ('" + tS|H (12) 

n (2b ZN { C F N f 7 C F N f 



r 2 = — f V L h 2 -h 3 + '—hi (13) 

6 V 8 4N C Nl 8 N c 12N C J V ; 

with hi = |j, h 2 = 67 gg 7 " 2 , h 3 = 47r2 24 ; 15 and h 4 = ^, also 3NLLO results have been 
derived [19 . Results from these approximations ^31 are shown in Fig. 1 together with the 
numerical solution of the MLLA evolution equations obtained in 1998 0] which takes into 
account all higher order corrections from this equation and fulfils the (non-perturbative) 
boundary condition JHJl. All curves are absolute predictions, as the parameter A (and 
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Qo in case of the numerical calculation) is adjusted from the growth of the total particle 
multiplicity in the e + e~ jets. The slow convergence of this ^fal expansion can be seen 
and there are still considerable effects beyond 3NLLO. The numerical solution is also in 
close agreement with the MC result at the parton level obtained [201 from the HERWIG 
MC above the jet energy Ej et > 15 GeV (£j e t = Q/2 in e + e~ annihilation) and ~ 20% 
larger at E- ]ct ~ 5 GeV. This overall agreement suggests that the effects not included in 
the MLLA evolution equation, such as large angle emission, are small. 

These numerical results are also compared in Fig. 1 with data from OPAL [201 where 
the data on gluon jets are derived from 3-jet events in e + e _ -annihilation. Note also that 
a proportionality constant relating partons and hadrons according to LPHD drops in 
the ratio r. Recent results from DELPHI ^] fall slightly below the curve by about 
20% at the lowest energies but converge for the higher ones; the CDF collaboration 
comparing quark and gluon jets at high pj> in pp collisions |22| finds the ratio r in the 
range 5 < E- ]ot < 15 GeV a bit larger, closer to the 3NLLO prediction, but with larger 
errors and therefore still consistent with the LEP results. 

An alternative calculation is based on the colour dipole model which treats the evo- 
lution of dipoles in NLL approximation and includes recoil effects |23| . It provides a 
good description of the data but includes an extra (non-perturbative) parameter which 
allows to adjust a low energy point. Whether such a non-perturbative input is definitely 
required by the DELPHI data depends in particular on the theoretical uncertainty in the 
extraction of N g from 3-jet events. 



5 Multiplicity moments of hadrons and sub-jets 

We now differentiate between sub-jets and hadrons in a jet of hard scale E = Q/2 in e + e~ 
annihilation. Sub-jets are defined by the resolution scale Q c [kr > Qc > A), hadrons are 
related to partons at scale Qo (kr > Qo) an d experience shows that Qo ~ A. 

In the DLA the multiplicity of partons at resolution Q c in a jet can be obtained 
analytically from J7J) with JSJ in terms of modified Bessel functions 

N g (Y) = p^Y + \ c {K {f3^fK)h{(3^Y + A c ) + I^^K^p^K)}. (14) 

where (3 2 = ^f^- with b = ^N c - |n/ and A c = ln^. At high energies I n ,K n rise 
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Figure 2: Partem multiplicity N in MLLA in a single gluon jet vs energy for fixed 
Qc = Qo (representing hadrons, full line) and at fixed energy Yq = \n(E/Qo) = 5 (LEP-1 
energy) for variable jet resolution Q c (l.h.s.) and the factorial moments F q vs. energy Y 
(r.h.s.). Numerical solutions of evolution eq., taken from 0. 

exponentially and one obtains in two limits for resolution Q c 

at high resolution (Q c -> Q ) : N ~ {[3 2 Y) 1/A In ( — ) exp y/(3(Y + A c ) (15) 

at low resolution (Q c — * E) : N — >• 1, (16) 

where we also used Kq(z) ~ ln(2/z) for small z. At high resolution for Q c — > A (A c — > 0) 
the parton multiplicity diverges logarithmically because of the Landau pole appearing 
in the running coupling. The pole is shielded by the cut-off Q c — Qq and at this value 
the parton multiplicity N reaches the hadron multiplicity Nh according to the LPHD 
prescription (up to an overall constant K). 

The multiplicity in MLLA comes out considerably smaller than in DLA, however the 
dependences on E and Q c are qualitatively the same. It is shown in Fig. [21 for fixed 
Qc = Qo the multiplicity N starts with N = 1 at threshold and rises ~ exp(cvln E); 
at fixed energy E we find N — 1 for Q c = E according to I jlfijl whereas N rises rapidly 
for Q c — > A and ultimately approaches the upper curve at Q c = Qo > A. The splitting 
between the upper and lower curve is entirely due to the running of the coupling a s . 

The full numerical solution of the MLLA evolution equations has been studied already 
in 0] . Similar results are obtained from the ARIADNE MC at the parton level with 
readjusted two parameters A, Qq according to the duality approach ( "APJADNE-D" ) . 
They are shown in Fig. in comparison with the experimental hadron multiplicities 
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Figure 3: Multiplicity N of hadrons (taken as N c h x 1.25) and multiplicity of jets in e + e 
annihilation at three energies Q together with parton Monte Carlo results (parameters 
A, Q fitted) as function of Y cs = \n(Q 2 /{Q 2 c + Q§)), Fig. from 0. 



(upper curve) and the jet multiplicities in e + e~ annihilation represented as superposition 
of 2 single jets (N = 2 at threshold) . The normalization of the hadron data is adjusted to 
the calculation; taking into account that the data do not include neutrals one finds that 
the proportionality factor K between multiplicities of hadrons and partons at scale Qq is 
close to K = 1 in agreement with previous findings 0] . In the comparison between data 
and calculations we used the variable Y cs ; it takes into account that in the experimental 
(and MC) jet algorithm the hadrons are resolved at Q c — but in the analytical MLLA 
calculation at Q c = Qq. 

A satisfactory overall description of the data can be obtained, especially in the tran- 
sition region between jets and hadrons, with parameters A = 400 MeV and Qq = 404 
MeV. Some disagreement with data of the jet curves occurs around Q c > 10 GeV which 
may be related to bb production, not included in the calculation. The differences between 
the various curves are determined by the behaviour of the running coupling: for fixed 
coupling all curves would coincide and behave like a power N cx (Q/Q c ) a - The rapid 
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Figure 4: Ratio of moments H q for hadrons as function of energy Q = 2E (fixed Qo) 
and for jets as function of resolution Q c at Q = 91 GeV using the variable Y cs as in Fig. 
m comparison with ARIADNE-D MC, from [5]. 

variation of the jet curves at their upper end comes from the closeness of the parameters 
Qo and A. 

Next we turn to the higher multiplicity moments. The energy dependence of the 
factorial moments F q in MLLA is also shown in Fig. |2 The calculation takes into ac- 
count energy conservation, therefore the threshold for production of q particles is shifted 
towards E t hr = qQo- Remarkably, all F q curves cross to good approximation at F q = 1, 
a Poissonian point. For a Poissonian, the kumulant moments K q and therefore the ratios 
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H q = K q /F q vanish for q > 1. At threshold one finds the rapid oscillations with q of 
the kumulants K q = (— l) q ^ 1 (q — 1)!, they continue up to the Poissonian point with 
decreasing amplitudes. At higher energies, the oscillation length of K q , H q increases and 
should approach asymptotically the DLA limit H q ~ 1/q 2 with all kumulants positive. 

In Fig. 01 we show as representative examples the moment ratios H2 and H3 for 
hadrons and jets as function of energy or resolution, respectively. One observes the 
approximate coincidence of the zeros of H2 and the minima near zero of H3 corresponding 
to the "Poissonian point" and the alternating signs below and same sign above this 
point. Again, the agreement of the structures in the region of highly resolved jets is well 
reproduced and the same is observed for all available higher moments jjj]. As function 
of order q one obtains for both hadrons and jets an oscillation pattern which depends 
on energy and resolution. The main characteristics can be derived from the numerical 
solution of the MLLA evolution equation and in good overall quantitative agreement 
with the ARIADNE-D MC. 



6 Conclusions 

The perturbative approach to multiparticle production based on the MLLA is found to 
work well, also for the correlation phenomena discussed here, both for hadrons and jets 
at variable resolution and it properly distinguishes between quark and gluon jets. An 
important condition for this success is the high accuracy of the calculation. The DLA 
at realistic energies does not always give qualitatively correct results and misses, for 
example, the Poissonian point and the oscillations of H q at higher energies. Also some 
additional terms in a a/oJ expansion are insufficient for a quantitative description of 
Ng/N q and the higher multiplicity moments. On the other hand, numerical solutions of 
the MLLA equation and the parton MC considered here come to a rather close description 
in terms of only two essential parameters A, Qq. 

The normalization parameter is found close to K = 1 4 which means that the global 
hadronic observables considered here can be described after replacing a hadron by a 
parton at resolution Qo > A. This description does not take into accont local effects like 
resonance production, but its success is in support of a dual description of a large class 
of hadronic and partonic observables. 
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